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ABSTRACT 

The apparently single WN 6 type star WR 134 (HD 191765) is distinguished among 
the Wolf-Rayet star population by its strong, presumably cyclical (V ~ 2.3 day), spectral 
variations. A true periodicity — which is still very much debated — would render WR 
134 a prime candidate for harboring either a collapsed companion or a rotating, large- 
scale, inhomogeneous outflow. 

We have carried out an intensive campaign of spectroscopic and photometric mon- 
itoring of WR 134 from 1989 to 1997 in an attempt to reveal the true nature of this 
object. This unprecedentedly large data set allows us to confirm unambiguously the 
existence of a coherent 2.25 ± 0.05 day periodicity in the line-profile changes of He II 
A4686, although the global pattern of variability is different from one epoch to another. 
This period is only marginally detected in the photometric data set. 

Assuming the 2.25 day periodic variability to be induced by orbital motion of a 
collapsed companion, we develop a simple model aiming to investigate (i) the effect of this 
strongly ionizing, accreting companion on the Wolf-Rayet wind structure, and (ii) the 
expected emergent X-ray luminosity. We argue that the predicted and observed X-ray 
fluxes can only be matched if the accretion on the collapsed star is significantly inhibited. 
Additionally, we performed simulations of line-profile variations caused by the orbital 
revolution of a localized, strongly ionized wind cavity surrounding the X-ray source. 
A reasonable fit is achieved between the observed and modeled phase-dependent line 

Dcpartcment de Physique, Universite de Montreal, CP. 6128, Succ. Centre- Ville, Montreal, Quebec, Canada, H3C 3J7; and Observa- 
toire du Mont Megantic; morel, sergey, moffat, eversber, stlouis@astro.umontreal.ca. 

* Present adress: Astrophysics Group, Imperial College of Science, Technology and Medicine, Blackett Laboratory, Prince Consort 
Road, London, SW7 2BZ, UK; email: morel@ic.ac.uk. 

** Present adress: Feinfocus Medizintechnik GmbH, Im Bahlbrink 11-13, 30827, Garbsen, Germany; email: t_eversberg@feinfocus.com. 
Dcpartamcnto de Astronomi'a, Universidad de Guanajuato, Apdo. Postal 144, 36000 Guanajuato Gto, Mexico; 
eenens@carina. astro, ugto. mx. 

Killam Research Fellow of the Canada Council for the Arts. 

Instituto de Astronomi'a, UNAM, Apdo. Postal 70-264, Mexico D.F. 04510, Mexico; gloria@astroscu.unam.mx. 
Sternberg Astronomical Institute, Universiteskij Prospect 13, 119899, Moscow, Russia; igor@sai.msu.su. 
Instituto de Astronomia, Apdo. Postal 877, C. P. 22860, Ensenada, B. C, Mexico; gag@bufadora.astrosen.unam.mx. 
McDonald Observatory, HET, P. O. Box 1337, Fort Davis, TX; grant@astro.as.utexas.edu. 

Instituto Nacional de Astrofi'sica, Optica y Electronica, Apdo. Postal 51, Puebla, Pue. 72000, Mexico; ocardona@inaoep.mx. 



profiles of He II A4686. However, the derived size of the photoionized zone substantially 
exceeds our expectations, given the observed low-level X-ray flux. 

Alternatively, we explore rotational modulation of a persistent, largely anisotropic 
outflow as the origin of the observed cyclical variability. Although qualitative, this 
hypothesis leads to greater consistency with the observations. 

Subject headings: stars: individual (WR 134) — stars: mass loss — stars: Wolf-Rayet 
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1. Introduction 

The recognition in the 1980s that some apparently 
single Wolf-Rayet (WR) stars exhibit seemingly peri- 
odic line profile and/or photometric variations argued 
for the existence of systems made up of a WR star 
and a collapsed companion (hereafter WR + c; where 
c stands either for a neutron star or a black hole), 
as predicted by the general theory of massive close- 
binary evolution (e.g., van den Heuvel & de Loore 
1973): O + O^WR+O^c + O^c + WR^c 
+ c. Due to the recoil of the first supernova explosion 
leading to c + O, the system generally acquires a high 
systemic velocity. If massive enough, the secondary 
evolves in turn into a WR star, at which point the 
system may have reached an unusually high galactic 
latitude for a population I star. These two peculiar- 
ities (runaway velocity and position), along with the 
existence of a surrounding ring nebula (presumably 
formed by matter ejected during the secondary mass 
exchange), were among the criteria initially used to 
select WR + c candidates (van den Heuvel 1976; Mof- 
fat 1982; Cherepashchuk & Aslanov 1984) .0 

A major breakthrough in the qualitative scenario 
described above comes from a recent redetermination 
of the distribution of the radio pulsar runaway ve- 
locities. In sharp contrast with earlier estimations 
(100-200 km s -1 ; Gunn & Ostriker 1970; Lyne, An- 
derson, & Salter 1982), these new values imply a mean 
pulsar kick velocity at birth of about 450 km s _1 
(Lyne & Lorimer 1994; Lorimer, Bailes, & Harrison 
1997). Although such a high kick velocity imparted 
at birth is still much debated (e.g., Hansen & Phinney 
1997; Hartman 1997), this result tends to show — as 
suggested by the apparent paucity of "runaway" OB 
stars with compact companions (Kumar, Kallman, & 
Thomas 1983; Gies & Bolton 1986; Philp et al. 1996; 
Sayer, Nice, & Kaspi 1996) — that the number of 
systems which would survive the first supernova ex- 
plosion is probably considerably lower than initially 
thought (see Brandt & Podsiadlowski 1995 vs Hellings 
& de Loore 1986). Indeed, the incorporation of up- 
to-date physics of supernova explosions in the most 
recent population synthesis models of massive bina- 

9 Since then, only Cygnus X-3 has been shown to probably belong 
to this class (van Kerkwijk et al. 1996; Schmutz, Geballe, & 
Schild 1996; but see Mitra 1998). Besides this system, the 
single-line WN 7 star WR 148 (HD 197406) constitutes one of 
the most promising candidates (Marchenko et al. 1996a). Two 
candidates for WR + c in 30 Doradus have also been proposed 
by Wang (1995). 



ries leads to a very small number of observable WR 
+ c systems in the Galaxy (N <^ 5; De Donder, Van- 
beveren, & van Bever 1997). This number contrasts 
significantly with the number of observationally se- 
lected candidates (N « 15; Vanbeveren 1991), espe- 
cially considering the necessarily incomplete nature 
of this sample. At face value, this discrepancy might 
imply that most of these WR + c candidates are spu- 
rious or, more interestingly, that other physical mech- 
anisms are at work to induce the large-scale, periodic 
variability inherent to some objects. 

At least in O stars, the progenitors of WR stars, 
it appears that large-scale periodic variability may be 
induced by rotating aspherical and structured winds 
(e.g., Fullerton et al. 1997; Kaper et al. 1997). The 
incidence of asymmetric outflows among the WR star 
population is probably lower, although recent spec- 
tropolarimetric (Harries, Hillier, & Howarth 1998), 
UV (St-Louis et al. 1995), optical (Marchenko et al. 
1998a), and radio observations (White & Becker 1996; 
Williams et al. 1997) have now independently re- 
vealed this peculiarity in a substantial number of WR 
stars. In order to account for the periodic variability 
inferred for the suspected WR + c candidates, the ro- 
tational modulation of a persistent, largely inhomo- 
geneous outflow could thus constitute in some cases 
an attractive alternative to the binary scenario, and 
would be consistent with the lack of strong, accretion- 
type X-ray emission (Wessolowski 1996). 

The nature of the variable WN 5 star EZ CMa 
(HD 50896; WR 6), which is generally considered 
as a prime candidate for an orbiting collapsed com- 
panion (Firmani et al. 1980), has been recently re- 
investigated by means of extensive spectropolarimet- 
ric (Harries et al. 1999), UV (St-Louis et al. 1995), 
and optical studies (Morel, St-Louis, & Marchenko 
1997; Morel et al. 1998). These studies reveal that the 
3.77 day period displayed by this object is more likely 
induced by the rotational modulation of a structured 
wind. Hints of long-term (from weeks to months) 
wind variability triggered by a spatial dependence of 
the physical conditions prevailing in the vicinity of the 
hydrostatic stellar "surface" were also found (Morel 
et al. 1998). 

Besides the search for WR + c systems, the mys- 
terious origin of such deep-seated variability (possi- 
bly induced by pulsational instabilities or large-scale 
magnetic structures in the case of EZ CMa) has 
prompted us to carry out intensive spectroscopic and 
photometric observations of other WR stars suspected 
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to exhibit such periodic line-profile variations (LPVs). 
In this context, WR 134 is regarded as a natural tar- 
get. 

2. Observational Background 

WR 134 is a relatively bright (v « 8.3) WN 6 star 
surrounded by a ring-type nebula embedded in the H 
II region S 109 (Crampton 1971; Esteban & Rosado 
1995). 

Following the discovery of spectacular line-profile 
and photometric changes intrinsic to this object (Bappu 
1951; Ross 1961), a major effort was directed to estab- 
lish the periodic nature of these variations. Cherepashchuk 
(1975) was the first to investigate photometric vari- 
ability of WR 134 on different timescales. Irregu- 
lar night-to-night light fluctuations (with rms am- 
plitude a «s 0.014 mag) were observed, without ev- 
idence for rapid (~ hourly) or long-term (~ monthly) 
changes. The first claim of periodic variability with V 
= 7.44 ± 0.10 day was made by Antokhin, Aslanov, 
& Cherepashchuk (1982) on the basis of a two-month 
interval of broadband photometric monitoring. The 
authors also reported small amplitude radial velocity 
variations present in simultaneously acquired spectro- 
scopic data (Kwr ~ 20-40 km s _1 ), consistent with 
the above period. This period was subsequently im- 
proved to V = 7.483 ± 0.004 day by Antokhin & 
Cherepashchuk (1984). Zhilyaev & Khalack (1996) 
established that the short-term stochastic variability 
presented by WR 134 is not related to this period, 
as it would be in the case of an orbiting collapsed 
companion. Moffat & Shara (1986) were unable to 
identify this period in their photometric data, al- 
though they found evidence for a 1.8 day periodicity, 
first tentatively reported by Lamontagne (1983) on 
the basis of radial velocity measurements. Following 
the recognition that the level of continuum flux from 
WR 134 is also irregularly variable on a timescale of 
weeks to months (Antokhin & Volkov 1987), periodic 
changes of the equivalent widths of He II A4859 with 
V = 1.74 ± 0.38 day were reported by Marchenko 
(1988). Then, Robert (1992) reported the existence 
of a 2.34 day periodicity in radial velocity variations, 
which could be a one-day alias of the 1.8 day pe- 
riod claimed in previous studies. A two-day quasi- 
periodicity in photometric data was also reported by 
Antokhin et al. (1992). Although not obvious in the 
spectroscopic data set obtained by Vreux et al. (1992) 
because of an unfortunate correspondence to a badly 



sampled frequency domain (see also Gosset, Vreux, 
& Andrillat 1994; Gosset & Vreux 1996), strong hints 
toward the reality of this two-day recurrence timescale 
come from the analysis of an extensive set of optical 
spectra by McCandliss et al. (1994), who derived a 
2.27 ± 0.04 day periodicity in the centroid, second 
moment and skewness of the line-profile time series. 
A relatively large-scale, likely cyclical, shift of excess- 
emission components superposed on the underlying 
line profiles was also observed. 

WR 134 also shows strong variations in broad- 
band polarimetry (Robert et al. 1989) and spectropo- 
larimetry (Schulte-Ladbeck et al. 1992). 

Numerous models have been put forward to ac- 
count for the intricate variability pattern, including: 
the presence of an orbiting neutron star (Antokhin & 
Cherepashchuk 1984), a disk connected to the cen- 
tral star by ever-changing filaments (Underhill et al. 
1990), or a bipolar magnetic outflow (Vreux et al. 
1992). 

No consensus has been reached yet concerning the 
possible periodic nature of the variations in WR 134. 
We present in this paper the analysis of a long-term 
campaign of optical photometric and spectroscopic 
monitoring, in an attempt to shed new light on this 
issue. Preliminary results concerning the photometric 
data subset were presented by Moffat & Marchenko 
(1993a, b) and Marchenko et al. (1996b). 

3. Observations and Reduction Procedure 
3.1. Spectroscopy 

WR 134 was observed spectroscopically between 
1992 and 1995. A journal of observations is presented 
in Table 1, which lists an epoch number, the dates 
of the spectroscopic observations, the interval of the 
observations in Heliocentric Julian Days, the obser- 
vatory name, the number of CCD spectra obtained, 
the selected spectral domain, the reciprocal dispersion 
of the spectra, and the typical signal-to-noise ratio 
(S/N) in the continuum. The spectra were reduced 
using the IRAF^| data reduction packages. Bias sub- 
traction, flat-field division, sky subtraction, extrac- 
tion of the spectra, and wavelength calibration were 
carried out in the usual way. Spectra of calibration 

'IRAF is distributed by the National Optical Astronomy Obser- 
vatories, operated by the Association of Universities for Re- 
search in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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lamps were taken immediately before, and after the 
stellar exposure. The rectification of the spectra was 
carried out by an appropriate choice of line-free re- 
gions, subsequently fitted by a low order Legendre 
polynomial. 

In order to minimize the spurious velocity shifts 
induced by an inevitably imperfect wavelength cali- 
bration, the spectra were coaligncd in velocity space 
by using the interstellar doublet Na I AA5890, 5896, 
or the diffuse interstellar band at 4501 A as fiducial 
marks. 

Each spectrum obtained during epochs I and II 
generally consists of the sum of three consecutive ex- 
posures typically separated by about 10 minutes, with 
no apparent short-term variability. This is expected 
in view of the longer timescales required to detect any 
significant motion of stochastic, small-scale emission- 
excess features travelling on top of the line profiles 
induced by outwardly moving wind inhomogeneities 
(Lcpinc & Moffat 1999). 

No attempts have been made to correct for the 
continuum level variability, owing to its relatively low 
amplitude and irregularity (§ 4.1.2). 

3.2. Photometry 

Previous work led us to expect a fairly complex 
light curve behavior. Thus, in attempting to reveal 
regular (i.e., periodic) variations, it was deemed es- 
sential to combine multi-epoch observations of high 
quality. We therefore organized a long-term cam- 
paign of UBV photometric monitoring of WR 134 
in 1990-1997 using the 0.25 m Automatic Photomet- 
ric Telescope (APT; Young et al. 1991). The APT 
data collected during the first 3 years of monitoring 
have been discussed by Moffat & Marchenko (1993a, 
b). The data generally consist of systematic, ~ day- 
to-day photometry with 1-2 measurements per night, 
each of accuracy of a w 0.005-0.008 mag. HD 192533 
and HD 192934 were used as check and comparison 
stars, respectively. 

In 1992 and 1995, we also organized additional 
multi-site, very intensive photometric campaigns, in 
an attempt to support the simultaneous spectroscopy. 
In 1992, we used the 0.84 m telescope of San Pe- 
dro Martir Observatory (Mexico) and the 0.6 m tele- 
scope of Maidanak Observatory (Uzbekistan, fSU), 
to observe WR 134 in narrow and broadband fil- 
ters. In Mexico, the star was observed with a two- 
channel (WR and guide star) photometer in rapid (5 



s resolution) or ultra-rapid (0.01 s; see Marchenko 
ct al. 1994) mode, through filters centered at A c 
= 5185 A (FWHM = 250 A, ~ stellar continuum), 
and A c = 4700 A (FWHM = 190 A, He II A4686 
emission line) . By subsequently re-binning the rapid- 
photometry data into 0.1 day bins, we achieved an ac- 
curacy of 0.002-0.003 mag per bin. Additionally, we 
used the two-channel device as a conventional one- 
channel photometer to obtain differential UBV ob- 
servations 1-2 times per night (individual accuracy of 
0.006-0.007 mag), with HD 192533 and HD 192934 
as a check-comparison pair. The single-channel pho- 
tometer of Maidanak Observatory was equipped with 
a single A c = 6012 A (FWHM = 87 A, ~ stellar con- 
tinuum) filter. HD 191917 was used as a comparison 
star, resulting in a mean accuracy of 0.003 mag per 
0.1 day binned data point. 

In 1995, we used the 0.84 m telescope of San Pe- 
dro Martir Observatory and one-channel photometer 
with a broadband V filter in differential photome- 
try mode, providing a typical accuracy of 0.003-0.005 
mag per 0.1 day binned point. WR 134 was also 
monitored by the Crimean Observational Station of 
GAISH (Crimea, Ukraine) with a one-channel pho- 
tometer (V filter) attached to a 0.6 m telescope, pro- 
viding a mean accuracy of 0.007 mag per 0.1 day bin. 

We also used all the photometric data secured by 
the HIPPARCOS astromctric satellite in 1989-1993 
(broadband Hp system; see Marchenko et al. 1998b). 
To match the accuracy of the ground-based observa- 
tions, we rebinned the HIPPARCOS data set to 1 
day bins, achieving a « 0.006 mag per combined data 
point. All relevant information concerning the pho- 
tometry is provided in Table 2. 

4. Results 

4.1. Search for Short-Term Periodicity 

Jf.A.1. Spectroscopy 

Inspired by the previous clear indication of a 2.3 
day periodicity in the line-profile variations found by 
McCandliss et al. (1994), we have performed a sim- 
ilar search in the present data set by calculating the 
power spectra (PS) using the technique of Scargle 
(1982) on the skewness, centroid, and FWHM time 
series of He II A4686, the strongest line. The cen- 
troid and skewness were calculated as the first mo- 
ment, and the ratio of the third and the (3/2 power 
of the) second central moments of the line profile, 
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respectively. In order to minimize the contributions 
of blends, both measurements were restricted to the 
portion of the profile above two in units of the con- 
tinuum. The FWHM was determined by a Gaussian 
fit to the entire line profile. A subsequent correction 
of the frequency spectrum by the CLEAN algorithm 
was performed in order to remove aliases and spurious 
features induced by the unevenly spaced nature of the 
data (see Roberts, Lehar, & Dreher 1987). The period 
search was performed up to the Nyquist frequency for 
unevenly spaced data: vn = (2At m j„,) _ . However, 
since no evidence for periodic signals was found at 
high frequencies, we will only display in the follow- 
ing the PS for frequencies up to vn — (2At mea n) ■ 
The data acquired during epoch III suffer from large 
time gaps and are inadequate in the search for peri- 
ods of the order of days. Also, because of the larger 
time span of epoch I compared to epoch II (and of 
the related higher frequency resolution), the period 
determination was mainly based on the former data 
set. The "raw" and CLEANed PS of the skewness, 
centroid, and FWHM time series of He II A4686 are 
presented in Figure 1. Note that the centroid vari- 
ations are mainly due to changes in the line profile 
morphology and thus do not reflect a global shift of 
the profile. The significance of the peaks in the PS 
can be estimated by means of the 99.0 % and 99.9 
% thresholds (Fig.l), giving the probability that a 
given peak is related to the presence of a determinis- 
tic signal in the time series (Scargle 1982). Note that 
these levels are only indicative in the case of unevenly 
spaced data (sec, e.g., Antokhin et al. 1995). 

The "raw" PS display numerous peaks, most of 
them being induced by the fairly complex temporal 
window and/or the presence of noise. Very few sig- 
nificant peaks remain in the CLEANed PS. In partic- 
ular, the highest peaks in the skewness and centroid 
PS are found at v « 0.442 and 0.447 d" 1 (with a typ- 
ical uncertainty of about 0.010 d _1 , as determined by 
the FWHM of the peaks), which translates into peri- 
ods of 2.26 ± 0.05 and 2.24 ± 0.05 day, respectively. 
These values are equal, within the uncertainties, to 
the period proposed by McCandliss et al. (1994): 
2.27 ± 0.04 day. The PS for the FWHM time se- 
ries presents a prominent peak at v\ rj 0.884 d _1 , 
which can very likely be identified with the first har- 
monic of the "fundamental" frequency vq suggested 
above (we adopt vq = 0.444 d _1 in the following). A 
highly significant signal is also found at v 2 ~ 0.326 
d" 1 (V = 3.07 ± 0.10 day); we will return to this 



point below. Note the absence of any trace of the 
7.483 day period (y w 0.134 d _1 ) proposed by An- 
tokhin & Cherepashchuk (1984). The periodic nature 
of the skewness, centroid, and FWHM variations with 
a frequency of v$ is confirmed — at least for epochs I 
and II — when the data are plotted as a function of 
phase (Fig. 2)]^] We have also calculated the equiva- 
lent widths (EWs) of He II A4686 by integrating the 
line flux in the interval 4646-4765 A. Although signif- 
icant, the EW variations show no clear evidence for 
phase-locked variability (Fig. 2). 

After demonstrating that the integrated properties 
of the line profile are periodic in nature, it is natu- 
ral to ask whether the same recurrence timescale is 
present in the detailed LP Vs. To this end, we per- 
formed pixel-to-pixcl CLEANing of the epoch I and 
II data sets. All spectra were rebinned in a simi- 
lar manner and were used as independent time se- 
quences, thus creating 1024 CLEANed PS along the 
wavelength direction. This procedure confirms the 
presence of vq, along with its first harmonic at v\ w 
0.88 d _1 over the entire line profile of He II A4686 
during epoch I, along with weak traces of this pe- 
riod in He II A4542, N V AA4604, 4620 and N III 
A4640 (Fig.3). As in the PS of the FWHM time se- 
ries of epoch I (Fig.l), a periodic signal at v 2 ~ 0.326 
d _1 is also present in the pixel-to-pixel CLEANed PS. 
The occurence of the same signal at v 2 in the FWHM 
and pixel-to-pixel CLEANed PS raises the possiblity 
that this feature is real. However, vq (or its harmon- 
ics) is generally (within the uncertainties) recovered 
in the epoch II PS, contrary to v 2 . More importantly, 
the gray-scale plot of He II A4686 does not show a 
coherent pattern of variability when folded with v 2 , 
whereas it does when folded with vq (see below). Also, 
the phase diagram of the FWHM data of epoch I is 
considerably noisier when folded with this frequency. 
Further observations are essential to indicate whether 
this period is genuinely spurious. Although the LPVs 
of He II A4686 for epoch II are undoubtedly periodic 
with a frequency of Uq (see below), only a weak signal 
appears in the pixel-to-pixel PS at this frequency on 
the top of He II A4542 and He II A4686. This may 
be related to the more complex pattern of variability 
compared to epoch I. 

The next analysis we have performed consisted in 

11 Because our period of 2.25 ± 0.05 day is indistinguishable from 
the one derived by McCandliss et al. (1994), we adopt their 
ephemeris in the following: HJD 2,447,015.753 + 2.27 E. 
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grouping the spectra for each epoch into 0.02 phase 
bins. The individual, binned He II A4686 line pro- 
files minus their corresponding unweighted means for 
each epoch are arranged as a function of phase in the 
upper panels of Figures 4a-4c. This line transition 
was chosen because its LPVs are quite representative 
of (but stronger and clearer than) that of the other 
He II spectral features (§ 4.3). In view of the long 
time span of the observations (13, 4, and 63 cycles 
for epochs I, II, and III, respectively), He II A4686 
presents a remarkably coherent phase-related pattern 
of variability. However, it is noteworthy that signifi- 
cant cycle-to-cycle differences in the line-profile mor- 
phology are often found, a fact which is reflected in 
the gray-scale plots by discontinuities in the pattern 
of variability between consecutive bins (especially for 
epoch III; Fig. 4c). Four main factors can induce this: 
(i) lack of long-term coherency in the pattern of vari- 
ability; (ii) artificial loss of coherency induced by the 
uncertainty in the period; (iii) presence of additional 
small-scale profile variations created by shocked (and 
possibly turbulent) material carried out by the global 
stellar outflow, as observed in other WR stars (Lepinc 
& Moffat 1999); and (iv) continuum flux variations for 
which no allowance has been made; we deem this last 
factor as a less important effect because of the small 
amplitude of the changes (see below) . Concerning the 
first two points, note that the smaller the time span of 
the observations (Table 1), the higher the coherency. 

4-1-2. Photometry 

For the analysis of the broadband photometry, we 
only use the ^-filter data as being less contaminated 
by emission lines (about 7 % of the total V flux), along 
with all available mcdiumband visual-filter photome- 
try, in total 734 observations binned to 0.1 day (1 day 
for the HIPPARCOS data). Note that for WR 134, 
the small-scale variations observed in U, B, and V are 
fairly well-correlated in general (Moffat & Marchenko 
1993b). 

Combining all 1989-1997 data and reducing them 
to the APT photometric system, we concentrate on 
a search for relatively short periods; a detailed in- 
vestigation of the temporal behavior of the "secular" 
components (i.e., V « 625 day and V w 40 day vari- 
ations) being presented in Marchenko et al. (1996b) 
and Marchenko & Moffat (1998a). 

We constructed a CLEANed PS of the whole pho- 
tometric data set, but failed to find any significant 
peaks in the range of interest, i.e., around the ex- 



pected 2.3 day component. This is possibly due to 
the lack of coherency of the period on long timescales. 
The long-term component with V w 625 day com- 
pletely dominates the PS, followed by a less evident 
power peak at V w 40 day. Pre-whitening from both 
long-period variations does not improve the situation 
in the high-frequency domain. 

In the analysis of the data subsets, we have en- 
countered some limited success. The "subsets" are 
naturally- imposed segments (e.g., summer monsoon 
and winter gaps at the APT, short-term campaigns 
of intense monitoring), with gaps between the seg- 
ments exceeding, or comparable to, the length of the 
segment filled by observations. We obtain 14 such seg- 
ments, covering the period HJD 2,447,859-2,450,636, 
with the number of observations per segment varying 
from 9 to 161. The most frequently appearing detail 
in the CLEANed PS is clustered around v = 1.3-1.4 
d _1 (in 5 out of 14 PS). Twice, we detect a fairly 
strong signal at v « 0.9 d" 1 : HJD 2,448,896-2,448,976 
(43 observations), and HJD 2,450,522-2,450,636 (82 
observations). Concerning the most abundant data 
set secured during the 1992 campaign (161 data 
points distributed over 135 days; practically all of 
them are shown on Fig. 5), the expected signal at vq 
= 0.440 ± 0.018 d _1 is enhanced only during a few 
cycles, between HJD 2,448,818-2,448,824 (Fig.5). We 
have used the data shown in the middle section of 
Figure 5 to construct a folded light curve with V = 
2.27 day (ephemeris from McCandliss et al. 1994). 
Despite the fact that the periodic signal found during 
the interval HJD 2,448,818-2,448,824 is identical to 
the one found in the independently and practically si- 
multaneously acquired spectroscopic data set — thus 
demonstrating that this detection is unlikely to be 
fortuitous — , the very noisy appearance of the folded 
light curve (Fig. 6) shows that the detection of the V = 
2.27 day periodicity is only marginal in the photomet- 
ric data set. Adopting v a w 0.44 d _1 as the principal 
frequency, we may interpret the v w 0.9 and 1.3-1.4 
d _1 components found in the other subsets as the first 
and second harmonics, respectively. 

Another intensive multi-site monitoring campaign 
in 1995 reveals the presence of the long-term V « 
40 day cyclic component (Marchenko et al. 1996b, 
and our Fig. 7), masking the relatively weaker v\ w 
0.88 d _1 frequency, along with the possible initiation 
of V ~ 7-8 day variations (perhaps as in Antokhin 
& Cherepashchuk 1984) at HJD 2,449,920-2,449,948. 
The latter phenomenon points to the transient char- 
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acter of the V w 7-8 day variations. 

4.2. Similarities of the Variations Across a 
Given Line Profile 

In an attempt to objectively and quantitatively ex- 
amine a potential relationship in the pattern of vari- 
ability displayed by different portions of a given line 
profile, we have calculated the Spearman rank-order 
correlation matrices (see, e.g., Johns & Basri 1995; 
Lago & Gameiro 1998) whose elements r(i,j) give 
the degree of correlation between the line intensity 
variations at any pixels i and j across the line profile 
(these matrices are symmetric and a perfect positive 
correlation is found along the main diagonal where i 
= j). These matrices for He II A4686 are shown for 
all epochs in Figure 8, in the form of contour plots, 
where the lowest contour indicates a significant posi- 
tive or negative correlation at the 99.9 % confidence 
level. 

An inspection of these matrices allows one to draw 
the following general conclusions: (i) There are vari- 
ous regions in the auto-correlation matrices present- 
ing a significant positive or negative correlation; (ii) If 
the variations were principally induced by the changes 
in the continuum flux level (§ 4.1.2), one would ob- 
serve a tendency for a positive correlation over a sub- 
stantial fraction of velocity space; this is clearly not 
observed; (iii) The main diagonal is not a straight 
line with a width corresponding to the velocity res- 
olution (at most 100 km s _1 in our case), as ex- 
pected if two contiguous velocity elements would vary 
in a completely uncorrelated fashion, but is consider- 
ably wider. This suggests that the same pattern of 
variability affects simultaneously a fairly large veloc- 
ity range of the profile, as independently found by 
Lcpinc, Moffat, & Hcnriksen (1996), who derived a 
relatively large mean line-of-sight velocity dispersion 
for the emission subpeaks travelling across the line 
profile of WR 134 (cr^ ss 440 km s _1 ) compared to 
other WR stars in their sample (<7£ « 100 km s _1 ). 

If we turn to a detailed epoch-to-epoch comparison 
of the matrices in Figure 8, we find significant dif- 
ferences (note that the phase coverage is sufficiently 
similar for all epochs for a direct comparison to be 
made). For example, a positive correlation is ob- 
served for epoch I at (- 500, + 2400) and (+ 200, + 
1800) km s _1 , whereas the variations in these velocity 
ranges are rather negatively correlated for epoch II. 
The paucity of contours for epoch III indicates that 
no significant relationship existed for this epoch be- 



tween the changes presented by different portions of 
the He II A4686 line profile. This leads to the con- 
clusion that the global pattern of variability is likely 
to differ notably in nature on a yearly timescale, al- 
though the LPVs are coherent over shorter timescales 
when phased with the 2.27 day period (Figs.4a-4c). 

4.3. Similarities Between the Variations of 
Different Lines 

The similar qualitative behavior of different lines 
has already been emphasized in the past (McCandliss 
1988; Marchenko 1988; Vreux et al. 1992; McCandliss 
et al. 1994). Here, we re-address this point by cal- 
culating the degree of correlation between the LPVs 
at different projected velocities (referred to the line 
laboratory rest wavelength) in two given line profiles. 
The spectra obtained during epoch III are well suited 
for this purpose because of their wide spectral cover- 
age. The correlation matrices of He II A5412 with He 
II A4542, He II A4686, He II A4859, and the doublet 
C IV A5806 are presented for this epoch in Figure 9. 
As revealed by the prevalence of contours along the 
main diagonal, the LPVs presented by the He II lines 
are generally well correlated (this is especially true for 
the relatively unblended lines He II A5412 and He II 
A4686) , emphasizing that they vary in a fairly similar 
fashion. This is remarkable as the presence of blends 
and/or noise tends to mask any potential correlation. 
Only the blue wings of He II A5412 and C IV A5806 
are positively correlated for epoch III. The changes 
affecting He II A4686 and He II A4542 are also well 
correlated for epochs I and II. 

5. Discussion 

5.1. The Duplicity of WR 134 Questioned 

In view of the long-suspected association of WR 
134 with an orbiting collapsed companion (Antokhin 
et al. 1982; Antokhin & Cherepashchuk 1984), we 
discuss below the implications of these observations 
with respect to this scenario. 

5.1.1. The Expected Emergent X-ray Luminosity 

An interesting issue is whether the predicted X-ray 
luminosity produced by the accretion process (after 
allowance for wind absorption) can be reconciled with 
the X-ray observations of WR 134. 

First, if we associate the 2.27 day periodicity with 
orbital motion, we can obtain an estimate of the loca- 
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tion of the compact object in the WR wind. Assum- 
ing for simplicity a circular orbit (i.e., that the circu- 
larization timescale is shorter than the evolutionary 
timescale; Tassoul 1990), a canonical mass for the sec- 
ondary as a neutron star, Mx = 1.4 M Q , and a mass 
for WR 134, M* = 11 M© (Hamann, Koesterke, & 
Wessolowski 1995), we obtain an orbital separation 
of a w 17 Rq, or a w 6 WR core radii (Hamann et al. 
1995)Q 

Second, the total X-ray luminosity (in ergs s _1 ) 
produced by Bondi-Hoyle accretion of a stellar wind 
onto a degenerate object can be expressed by the fol- 
lowing relation (Stevens & Willis 1988): 



L 
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2.03 x 10 57 T] M M x a~ 2 v(a) 
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L E 



(1) 

where rj is the efficiency of the conversion of gravi- 
tational energy into X-ray emission; M the mass-loss 
rate of the primary (in M© yr _1 ); v(a) the WR wind 
velocity at the secondary's location (in km s _1 ), and 
Le the Eddington luminosity (in ergs s _1 ) given by 
Le ~ 5.02 x 10 37 Mx/cr e (the electron scattering co- 
efficient <r e = 0.35). In this expression, we neglect the 
orbital velocity of the secondary relative to the (much 
larger) WR wind velocity. For the latter, we start by 
arbitrarily adopting the well-established mean veloc- 
ity law for O stars of the form (Pauldrach, Puis, & 
Kudritzki 1986)F|: 



v(r) 



(1 - R*/r) 



0.8 



(2) 



Note that at the assumed location of the secondary 
(r w 6 R+; see above) the wind has almost reached its 
terminal velocity, Voo = 1900 km s _1 (Rochowicz & 
Niedzielski 1995). We further assume: M = 8 x 10~ 5 
M© yr^ 1 (Hogg 1989), and t] = 0.1 (McCray 1977). 
This yields: Lx ~ 1.4 x 10 37 ergs s _1 . 

We have modeled the emergent X-ray flux by a 
power law with energy index a, substituted above a 
characteristic high-energy cutoff E c by the function 
exp[(_E c — E)/Ef], as commonly observed in accre- 
tion powered pulsars (White, Swank, & Holt 1983; 



12 Since the observationally determined masses of WN 6 stars in 
binary systems show a large scatter (M* « 14 Mq in WR 153; 
St-Louis et al. 1988 — M* « 48 M in WR 47; Moffat et al. 
1990), we use here a value given by atmospheric models of WR 
stars (adopting M+ as 14 Mq or 48 Mq does not qualitatively 
change the conclusions presented in the following). 

13 We explore below the impact on Lx of more recently proposed 
(but not yet well established) revisions for the velocity law and 
mass-loss rate of WR stars. 



Krctschmar et al. 1997). The values of a, E c , and 
Ef are chosen to be roughly representative of accret- 
ing pulsars: a « - 0.2; E c w 15 keV; E f « 15 keV 
(White et al. 1983). The X-ray spectrum was normal- 
ized to the total X-ray luminosity determined above, 
i.e., L x = 1.4 x 10 37 ergs s -1 . 

The next step is to evaluate the attenuation of the 
beam of photons as they propagate through the stel- 
lar wind. Since the absorption properties of ionized 
plasmas are dramatically dependent on their ioniza- 
tion state (e.g., Woo et al. 1995), one has first to de- 
termine to what degree the accreting compact object 
ionizes the surrounding stellar wind. The presence of 
an immersed, strongly ionizing companion is expected 
to create an extended X-ray photoionized zone in its 
vicinity (Hatchett & McCray 1977). An investigation 
of the effect of the ionizing X-ray flux coming from 
the secondary on the surrounding wind material re- 
quires a detailed calculation of the radiative transfer 
(Kallman & McCray 1982). However, good insight 
into the ionization state of an optically thin gas illu- 
minated by a X-ray point source can be obtained by 
considering the quantity (Hatchett & McCray 1977): 
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with: n(r) the local number density of the gas; rx 
the distance from the X-ray source. This expression 
is derived using the mass continuity equation: M = 4 



7r r 2 m n(r) v(r). Here, m is the average mass per ion 
(we assume a pure helium atmosphere). The velocity 
of the material as a function of r is given by equa- 
tion (2). The variations of log£(r, r x ) are illustrated 
as viewed from above the orbital plane in Figure 10. 
In the limit £ — > ergs cm s _1 , the ionization bal- 
ance of the material is unaffected by the presence of 
the X-ray emitter, i.e., is principally governed by the 
radiation field of the WR star. We stress that this 
model assumes an optically thin plasma. In reality, 
the mean path length of the X-ray photons is consid- 
erably shorter, and it is thus likely that the X-ray pho- 
toionized zone would be dramatically less extended — 
especially in the direction toward the primary — than 
sketched here. 

The predicted (after wind absorption) X-ray lu- 
minosities for the wavebands corresponding to the 
PS PC and IPC detectors onboard the ROSAT and 
Einstein observatories (0.2-2.4 and 0.2-4.0 keV, re- 
spectively) have been computed for three illustrative 
cases spanning a wide range in the wind ionization 
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state, namely log£ = 2.1, 1.8, and ergs cm s -1 . The 
photoelectric atomic cross sections for ionized plas- 
mas are taken from Woo et al. (1995). Thompson 
scattering by the free electrons was also taken into 
account. The two former values, log£ — 2.1 and 1.8 
ergs cm s _1 , are roughly what would be expected in 
the framework of our model (Fig. 10), whereas in the 
latter case, log£ = ergs cm s -1 , we explore the 
(somewhat unrealistic) case in which the presence of 
the neutron star has no effect on the surrounding stel- 
lar wind. In this case, photoelectric cross sections for 
cold interstellar gas were used (Morrison & McCam- 
mon 1983). Because the WR wind material is, by 
nature, far from neutral, the derived X-ray luminosi- 
ties are here strictly lower limits. 

Pollock (1987) reported on two pointed IPC obser- 
vations of WR 134 obtained with the Einstein satel- 
lite, with no evidence for variability (see his Table 5). 
The derived (corrected for interstellar extinction) X- 
ray luminosity in the 0.2-4.0 keV band, 4.6 i 1.6 x 
10 32 ergs s _1 , being consistent with the upper limit 
derived by Sanders et al. (1985): 4.5 x 10 32 ergs 
s _1 (when scaled to the same adopted distance of 2.1 
kpc; van der Hucht et al. 1988). On the other hand, 
a single measurement (due to the ROSAT satellite) 
is available in the 0.2-2.4 keV range (Pollock, Haberl, 
& Corcoran 1995). With no clear indication of X- 
ray variability, we will assume in the following that 
the typical observed (and corrected for interstellar ex- 
tinction) X-ray luminosities from WR 134 are: 0.46 
± 0.22 and 4.6 ± 1.6 x 10 32 ergs s" 1 in the 0.2-2.4 
and 0.2-4.0 keV bands, respectively.^ 

These values can be directly compared with the 
predicted X-ray luminosities after wind absorption 
shown as a function of orbital phase in Figure 11. For 
log£ = 2.1 and 1.8 ergs cm s _1 , the observed luminosi- 
ties are systematically much lower than expected in 
the framework of our model, with a deficiency reach- 
ing 2-3 orders of magnitude. For log£ = ergs cm s _1 , 
the ROSAT data are consistent with the expectations. 
However, and unless the Einstein satellite observed 
WR 134 near X-ray eclipse, the same conclusion does 
not hold for the data in the 0.2-4.0 keV band, with 
an order of magnitude deficiency. Although mildly 
significant, this deficiency shows that, even in the ex- 
treme case in which the presence of the accreting neu- 

Unfortunately, the relatively large uncertainty in the period 
(§4.1.1) does not allow to fold in phase the ROSAT and Ein- 
stein data, thus preventing the construction of an observed 
X-ray light curve of WR 134. 



tron star has a negligible effect on the WR wind, the 
observed and predicted X-ray luminosities can hardly 
be reconciled. Overall, and although we stress that 
more detailed and rigorous calculations are necessary, 
the deficiency of the observed X-ray flux must thus be 
regarded as serious. This conclusion is bolstered if one 
considers that wind radiative instabilities may largely 
account for the observed fluxes (typically ifj 32 ~ 33 ergs 
s" 1 in the 0.2-2.4 keV range; Wessolowski 1996). 

Assuming a lower mass-loss rate due to clumping 
(Moffat & Robert 1994; Nugis, Crowther, & Willis 
1998) may, at best, reduce the discrepancy by one 
order of magnitude (see equation [1]). In contrast, 
however, assuming a "softer" v(r) law (i.e., (3 > 0.8; 
see equation [2]) for WR stars (Schmutz 1997; Lepine 
& Moffat 1999; Antokhin, Cherepashchuk, & Yagola 
1998), or taking into account the existence of a ve- 
locity plateau inside the photoionized zone (Blondin 
1994), exacerbates the discrepancy. 

Considering that the soft X-ray flux of WR 134 
is very similar to that of other bona fide single WN 
stars (Pollock 1987), and was even among the low- 
est found for single WN 6 stars during the ROSAT 
all-sky survey (Pollock et al. 1995), this deficiency 
of predicted X-ray flux suggests centrifugal (or mag- 
netic) inhibition to prevent the gravitational capture 
by the compact object of the primary's wind material 
(the so-called "propeller effect"; Lipunov 1982; Cam- 
pana 1997; Cui 1997; Zhang, Yu, & Zhang 1998). 
Note that orders of magnitude deficiency in the X- 
ray flux is observed in some high-mass X-ray binaries 
(HMXRBs; e.g., Taylor et al. 1996). 

5.1.2. Modeling the Line-Profile Variations Caused 
by an Ionizing X-ray Source 

Keeping in mind that the accretion process would 
be, as suggested above, much less efficient than ex- 
pected in the Bondi-Hoyle approximation — and 
therefore that the influence of the X-ray source on 
the surrounding wind material would be far less se- 
vere than sketched in Figure 10 — we performed sim- 
ulations of the LPVs caused by the orbital revolution 
of a localized, strongly ionized wind cavity. 

We proceed with numerical simulations of the ob- 
served LPVs of the representative He II A4686 line, 
binning the epoch I spectra to 0.1 phase resolution. 
We modify the SEI code (Sobolev method with ex- 
act integration: Lamers, Cerruti-Sola, & Perinotto 
1987; hereafter LCP) to allow for the variation of the 
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source function in a 3-D space, thus breaking the ini- 
tially assumed spherical symmetry of the WR wind. 
The variation of the unperturbed source function as 
a function of the radial distance from the star is de- 
scribed by equation (4) of LCP. The wind dynamics is 
assumed to be identical inside, and outside the pho- 
toionized cavity. Specifically, we adopt a /3-velocity 
law with an exponent (3 = 3 and a ratio of the wind 
velocity at the inner boundary of the optically thin 
part of the WR wind to the terminal velocity, wq, of 
0.25 (see equation [35] of LCP). 

We start from the assumption that the underly- 
ing, unperturbed line profile is not affected by any 
phase-dependent variations. We fit this reference pro- 
file with the standard SEI code, preserving complete 
spherical symmetry. Throughout these simulations, 
we adopt tt = 0.60, ot\ — 6.0, c*2 = 0.5, e' = 12.0, 
Bo = 1.2, <it = 0.2, and wd = 0.25 (these values 
are uncertain to about 30 %). We refer the reader 
to LCP for a complete description of these parame- 
ters. Note that these values slightly differ from those 
adopted by Marchenko & Moffat (1998b) due to a dif- 
ferent approach in modelling (adopting their values, 
however, would not qualitatively modify the conclu- 
sions drawn in the following). Note that while the 
issue of the flattened, asymmetric wind in WR 134 
is appealing (Schultc-Ladbeck et al. 1992; Moffat & 
Marchenko 1993b), this inevitably introduces a great 
and probably unnecessary complexity into the simu- 
lations. 

There are two obvious choices for the reference pro- 
file of the unperturbed wind: a minimum-emission or 
maximum-emission profile, derived from a smoothed 
minimum/maximum emissivity at a given wavelength 
for a given phase (Fig. 12). By introducing 3-D vari- 
ations into the source function of the unperturbed 
wind corresponding to the minimum-emission refer- 
ence profile, we immediately find that we are not 
able to reproduce the observed LPVs with any rea- 
sonable choice of the parameters. Thus we adopt the 
maximum-emission representation of the unperturbed 
profile as a starting approximation. 

To simulate the X-ray induced cavity in the other- 
wise spherically symmetric wind, we require the fol- 
lowing free parameters: — the azimuthal extension 
of the cavity; Az — Ap — line-of-sight and impact- 
parameter extension of the cavity; Azq — line-of-sight 
distance of the cavity's inner boundary from the stel- 
lar core; i — orbital inclination; k — coefficient de- 
scribing the deviation of the source function within 



the cavity, namely S (cavity) = k x S(wind); (f>o - 
cross-over phase (frontal passage of the cavity). We 
assume a circular orbit. A sketch of the adopted ge- 
ometry is shown in Figure 13. 

With these parameters, we attempt to fit all the 
phase-binned profiles, minimizing the observed minus 
model deviations in the x 2 sense, while reproducing 
the observed TVS of He II A4686 (Fig.12). We also 
attempt to make the (O - C) profile deviations dis- 
tributed as evenly as possible over the entire phase 
interval. We concentrate on the goodness of the fit 
for the velocities not exceeding ± 0.8 v^, since it is 
impossible to match the red wing of the He II A4686 
profile (electron scattering effect; see Hillier 1991) as 
well as the bluemost portion (partial blending with N 
III transitions). 

After an extensive search for optimal parameters, 
we find that the following set reproduces the observed 
LPVs reasonably well (Fig. 14): = 140 ± 10°, Az 
= Ap = 8 ± 0.5 i?*, Az = 2.7 ± 0.2 R+, i = 65 ± 
5°, k = 0.07 ± 0.01 (i.e., practically zero emissivity 
from the X-ray eroded wind), and <j>o = 0.55 ± 0.03. 
This value of (f>o means that the X-ray source passes 
in front when the continuum flux undergoes a shallow 
minimum (Fig. 6). The largest (O - C) deviations are 
concentrated around phases 0.0 and 0.5. By no means 
are we able to reproduce the extended emissivity ex- 
cesses at v w - 0.5 Woo [4> — 0.35) and v = + (0.5-0.7) 
^oo {4> = 0.65 and 0.75). These deviations might be 
generated by a low-amplitude velocity shift of the un- 
derlying profile as a whole, due to binary motion (par- 
tially inducing the centroid variations seen in Fig. 2). 
We will not consider this possibility further, as it will 
add more free parameters into the model. We are 
able to reproduce only the general appearance of the 
TVS as a structured, non-monotonic function with 3 
distinct maxima, without precise match of their am- 
plitudes and positions (Fig.12). 

According to Model 5 of Kallman & McCray (1982), 
helium is fully ionized for log £ ;> 1.5 ergs cm s" 1 . The 
derived azimuthal and radial extensions of the cavity 
indeed match fairly well the size of the highly ionized 
zone with log£ w 1.5 ergs cm s^ 1 (Fig. 10), especially 
if we take into account its large azimuthal extension. 
However, since the accretion is likely to be partially 
inhibited and the WR wind is not optically thin, the 
real photoionized cavity will be much less extended 
than sketched in Figure 10. For example, reduction of 
the efficiency of the X-ray generation by 1-2 orders of 
magnitude would dramatically reduce the dimension 
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of the log£ « 1.5 cavity (to the size of the log£ ;> 2.5 
cavities sketched in Figure 10; see equation [3]). The 
local nature of this putative zone is also suggested by 
an inspection of the P Cygni absorption profile vari- 
ations of He I A4471 and N V A4604. The epoch I 
spectra have been groupped into broad phase bins of 
width 0.2 (the P Cygni absorptions are very weak, 
thus requiring an extremely high S/N), and are over- 
plotted in Figure 15. It is apparent that potential 
phase-related variations do not significantly exceed 
the observationally- inflicted accuracy. Thus, neither 
the base of the wind (N V A4604), nor the relatively 
distant regions (He I A4471) arc seriously affected by 
the presence of the photoionized cavity. An inspec- 
tion of 16 archive IVE SWP spectra of WR 134 ob- 
tained during the period 1989 November 30 - Decem- 
ber 6 (see Table 3) also shows that no UV lines with 
well-developed P Cygni absorption components show 
any significant variability in their absorption troughs. 
Hints of relatively weak, possibly phase-locked varia- 
tions are only evident in the emission parts. 

We conclude that despite the encouraging similar- 
ity between the observed and modeled line-profiles, 
the relatively large derived size of the cavity (R « 
8 R+) along with its large azimuthal extension seem 
to enter in conflict with the low level of X-ray flux 
observed from WR 134. 

5.1.3. On the Origin of the Long-Term Modulations 

Another question we may ask is whether this model 
is able to account for the V w 40 and V w 625 
day recurrence timescales present in the photometric 
data (Marchcnko ct al. 1996b; Marchenko & Moffat 
1998a), as well as the epoch-dependent nature of the 
LPVs (Figs.2 and 4). 

As the wind approaches the compact object, it 
is gravitationally focused, forming a standing bow 
shock. This leads to the formation of an extended 
wake downstream (Blondin ct al. 1990). This process 
is suspected to be non-stationary in HMXRBs, with 
the possible formation of a quasi-cyclical "flip-flop" 
instability (Benensohn, Lamb, & Taam 1997; but see 
Ruffcrt 1997), with characteristic timescales compa- 
rable to the flow time across the accretion zone (« 
minutes); too short a timescale to account for any 
observed long-term variations in WR 134. 

Alternatively, it is tempting to associate the 40 day 
recurrence timescale to the precession of a tilted ac- 
cretion disk as suggested in some HMXRBs (Hcermskcrk 



& van Paradijs 1989; Wijnands, Kuulkers, & Smale 
1996). The condition for the creation of a persistent 
accretion disk around a wind-fed neutron star orbiting 
an early-type companion can be written as (Shapiro 
& Lightman 1976): 

2.7xlO- n ^B- 2 4/7 R- w/7 M™ /7 M- 4 p- 2 RiLV 7 > 1 

(4) 

With B12 the surface magnetic field of the neutron 
star in units of 10 12 Gauss, Rio the radius of the neu- 
tron star in units of 10 km, M x and M* in units of 
solar masses, P in units of days, i?* in units of solar 
radii, and Lx in units of ergs s _1 . 7 is expressed by 
the relation: 7 w C, 1 / 2 [v esc /v(a)j where £ is a dimen- 
sionlcss quantity accounting for the deviation from 
the Hoyle-Lyttleton treatment, and v esc the escape 
velocity at the WR "photosphere" (we neglect the or- 
bital velocity of the neutron star relative to the WR 
wind velocity). Taking B12 = Rio — ( = 1, wc obtain 
a value much below unity (~ 10~ 4 ), suggesting that 
accretion mediated by a persistent disk in WR 134 
would be very unlikely. 

A viable mechanism that might, however, account 
for the monthly changes in the spectroscopic pattern 
of variability (Figs.2 and 4) is to consider long-term 
changes in the accretion rate (thus varying the size of 
the ionized cavity surrounding the compact compan- 
ion). 

5.2. Wind-Related Variability? 

Although the encouraging agreement between the 
observed and modeled LPVs (Fig. 14) leads to the con- 
clusion that the model involving an accreting com- 
pact companion cannot be completely ruled out by 
the present observations, the very low observed X- 
ray flux of WR 134 which imposes inhibited accre- 
tion, as well as the long-term changes in the accre- 
tion rate required to account for the epoch-dependent 
nature of the LPVs, set serious constraints on this 
model. This, together with the inconclusive search 
for rapid periodic photometric variations, possibly in- 
duced by a spinning neutron star (Antokhin et al. 
1982; Marchenko et al. 1994) can be used to argue 
that the existence of such a low-mass collapsed com- 
panion is questionable. 

Alternatively, it is conceivable that we observe 
rotationally-modulated variability in WR 134. In- 
deed, the global pattern of spectral variability resem- 
bles in some aspects that of the apparently single WR 
star EZ CMa, for which the existence of a rotation- 
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modulated, structured wind has been proposed (St- 
Louis ct al. 1995; Morel et al. 1997, 1998; Harries 
et al. 1999). By analogy with EZ CMa, the He II 
line profiles vary in a fairly similar fashion (Vreux et 
al. 1992, and our Fig. 9), changes in line skewness 
or FWHM generally demonstrate a rather simple and 
well-defined behavior when phased with the 2.27 day 
period (Fig. 2), and a positive/negative correlation be- 
tween the pattern of variability presented by different 
parts of the same line profile is occasionally found 
(Fig. 8). Other outstanding properties shared by WR 
134 and EZ CMa are the epoch-dependent nature of 
the variations, as well as the substantial depolariza- 
tion of the emission lines (Harries et al. 1998); this 
last point being generally attributed to an equatorial 
density enhancement (see also Ignace et al. 1998). 

Another piece of evidence pointing to the similar- 
ity between WR 134 and EZ CMa comes from the 
broadband (~ continuum) polarimetric observations 
of Robert et al. (1989). Plots of the Stokes parame- 
ters Q and U of WR 134 versus phase are shown in 
Figure 16 (because of a likely loss of coherency over 
long timescales, the two data subsets separated by a 
year are plotted separately). The data show strong 
epoch-dependent variations, with a clear single-wave 
variation in U (less clear in Q) in 1985, and a fairly 
clear double-wave modulation in Q (less clear in U) 
in 1986. This behavior is qualitatively very similar to 
what is observed in EZ CMa (e.g., Robert et al. 1992), 
although the polarimetric variations are not as clearly 
phase-locked in WR 134. Such a drastic change from 
a single-wave to a double-wave modulation in data 
taken one year apart is not easily acommodated by 
a binary hypothesis. However, neither the amount, 
nor the time coverage of the polarimetric data allow 
to regard it as a decisive evidence favoring the single- 
star interpretation. Although the 7 nights of linear 
spectropolarimetry of WR 134 by Schulte-Ladbeck et 
al. (1992) do show variations, both in the continuum 
and in the lines, the data are sparsely spread out over 
6 months, making it impossible to separate variations 
on the 2.3-day cycle from long-term epoch-dependent 
variations. Certainly, further more intense spectropo- 
larimetry on different timescales in all 4 Stokes' pa- 
rameters will prove quite interesting. 

Relying on the EZ CMa and WR 134 similarities, 
one may speculate that the 2.3 day periodicity in 
WR 134 is induced by a small number of spatially 
extended, relatively long-living and rotating wind 
streams whose formation is possibly triggered by pho- 



tospheric perturbations, such as magnetic structures 
or pulsations (e.g., Cranmer & Owocki 1996; Kaper 
ct al. 1997). Contrary to EZ CMa, however, WR 134 
docs not show any detectable phase-related changes 
in the UV or optical P Cygni absorption troughs, nor 
a clear correlation between changes occuring near the 
hydrostatic stellar "surface" and in the wind. 

By speculating about the existence of (non)radial 
pulsations in WR 134 (although a search for related 
rapid light variations has been hitherto inconclusive; 
Chcrepashchuk 1975; Antokhin et al. 1992), and ex- 
ploring the idea of interactions between different pul- 
sational modes, one is led to the notion of quasi- 
periodic ejection of shell-like structures. This idea 
can be explored via the application of the model in- 
troduced above for the anisotropic wind (§ 5.1.2). We 
allow for the formation of a shell with enhanced op- 
tical depth in the otherwise spherically symmetric 
wind. The shell geometry may be variable but, for 
the sake of simplicity, we mainly explore the spher- 
ically symmetric case. The shell can slowly propa- 
gate outwards at a given rate, slightly expanding and 
gradually approaching the physical conditions in the 
unperturbed surrounding wind. Indeed, at least the 
variations in the blue wing of He II A4686 (Figs. 4a 
and 14) somewhat mimic the expected behavior of 
such an outwardly-moving structure. However, there 
is also pronounced asymmetry between the blue and 
red parts of the He II A4686 profile, especially obvious 
around <fi — 0.35-0.55 (Fig. 14). This asymmetry can- 
not be reproduced by any adjustment of the model 
free parameters, unless one makes the contrived as- 
sumption that the front- and back-side lobes of the 
shell are formed under different physical conditions. 
This leads us to discard the model of an expanding 
shell with central symmetry. 

Although the existence of a globally inhomogc- 
neous outflow in WR 134 is attractive, such an inter- 
pretation is also not without difficulties. For instance, 
the origin of the long-term (quasi)periodic changes in 
continuum flux, and of the epoch-dependency of the 
spectral changes remains unexplained. Magnetic ac- 
tivity, by inducing long-term changes in the global 
wind structure, may constitute a convenient, although 
largely ad hoc, way to accomodate this aspect of the 
variability. 
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TABLE 1 - JOURNAL OF SPECTROSCOPIC 
OBSERVATIONS. 







HJD 




Number 


Spectral 


Reciprocal Dispersion 




Epoch 


Date 


- 2440000 


Observatory" 


of Spectra 


Coverage (A) 


(A pix- 1 ) 


S/N 


I 


92 Jul - Aug 


8813-8843 


SPM 


65 


4380-4790 


0.41 


~ 160 


II 


93 Oct 


9258-9269 


SPM 


38 


3990-4840 


0.83 


~ 115 


III 


95 May - Oct 


9860-10005 


OMM, DAO 


32 


4360-5080 


1.62 


~ 210 










30 


5050-5945 


1.62 


~ 290 



a SPM: San Pedro Martir Observatory 2.1 m (Mexico); OMM: 
Obscrvatoire du Mont Mcgantic 1.6 m (Canada); DAO: Do- 
minion Astrophysical Observatory 1.2 m (Canada). 



TABLE 2 - JOURNAL OF PHOTOMETRIC 
OBSERVATIONS. 





HJD 




Number of Binned 


Mode of 


Filter 


Date 


- 2440000 


Observatory" 


Observations 


Observation 


Set 6 


89 Nov - 93 Feb 


7859-9045 


HIPPARCOS 


44 




Hp (~ B + V) c 


90 May - 97 Jul 


8025-10636 


APT 


483 


Diff. phot. 


V 


92 Jun - Jul 


8800-8828 


SPM 


28 


2-ch. phot. 


5185 (250) 


92 Jun - Jul 


8800-8832 


SPM 


23 


Diff. phot. 


V 


92 Jul 


8819-8833 


MO 


55 


Diff. phot. 


6012 (87) 


95 Jun - Jul 


9885-9903 


CR 


13 


Diff. phot. 


V 


95 Jun - Aug 


9891-9949 


SPM 


88 


Diff. phot. 


V 



a APT: Automatic Photometric Telescope 0.25 m (USA); SPM: 
San Pedro Martir Observatory 0.84 m (Mexico); MO: Maid- 
anak Observatory 0.6 m (Uzbekistan); CR: Observational Sta- 
tion in Crimea 0.6 m (Ukraine). 

b For the midband filters, we provide the central wavelength 

and FWHM (in A). 

c See van Leeuwen et al. (1997). 
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TABLE 3 - IUE SWP HIGH RESOLUTION SPECTRA OF 
WR 134. 



SWP Image Number 


JD - 2440000" 


Phase" 


37705 


7861.092 


0.396 


37707 


7861.213 


0.449 


37718 


7863.008 


0.240 


37725 


7863.981 


0.669 


37727 


7864.097 


0.720 


37734 


7865.136 


0.178 


37735 


7865.184 


0.199 


37737 


7865.315 


0.256 


37739 


7865.427 


0.306 


37741 


7865.543 


0.357 


37743 


7865.655 


0.406 


37745 


7865.765 


0.455 


37748 


7865.972 


0.546 


37750 


7866.095 


0.600 


37755 


7866.976 


0.988 


37757 


7867.103 


0.044 



a Julian Date at the midpoint 
of the exposure (2400 s for all spectra). 

b According to the ephemeris of 
McCandliss ct al. (1994): JD 2447015.753 + 2.27 E. 
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Fig. 1. — Power spectra (normalized to the total 
variance of the data) of the skewness, centroid, and 
FWHM time series of He II A4686 for epoch I. Dotted 
line: raw PS; solid line: CLEANed PS. The short- 
and long-dashed horizontal lines indicate the 99.0 % 
and 99.9 % confidence levels for the presence of a de- 
terministic signal in the time series, respectively. The 
number of iterations of the CLEAN algorithm, N, has 
been set to 3400 (skewness), 380 (centroid), and 4500 
(FWHM); a gain g = 0.2 has been used throughout. 
The algorithm is generally insensitive to the combi- 
nation of these two parameters in the range N ;> 100 
and g < 0.5-0.8. 
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Fig. 2. — Skewness and centroid variations, devia- 
tions of the FWHM around the mean value (in km 
s -1 ), and EW variations (normalized by division to 
the mean value) of He II A4686 for each epoch, as 
a function of phase. The 2-cr error bars for the 
EW values were calculated according to Chalabaev 
& Maillard (1983). As everywhere in this paper, 
the ephemeris of McCandliss et al. (1994) has been 
adopted: HJD 2,447,015.753 + 2.27 E. 
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4500 4600 

Wavelength (A) 



Fig. 3.— The pixel-to-pixel CLEANcd PS of the rec- 
tified spectra of epoch I. The lower and right-hand 
panels show the mean spectrum and the normalized, 
integrated power over all pixels, respectively. The 
highest peak corresponds to a frequency: v$ w 0.444 
d _1 . The gain and the number of iterations have been 
set to 0.3 and 250, respectively. 
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Fig. 4. — Gray-scale plots of the time series of the 
residuals of He II A4686 for each epoch, as a func- 
tion of phase. The residuals are the binned (to 0.02 
phase resolution) He II A4686 line profiles minus their 
corresponding unweighted means for each epoch. Ex- 
cess emission components appear brighter in these 
plots. The middle portion of each panel shows a su- 
perposition of the rectified profiles. The values of the 
temporal variance spectrum (TVS; Fullerton, Gies, 
& Bolton 1996), along with the horizontal dashed 
line indicating the 99.0 % confidence level for signifi- 
cant variability, are displayed in the lower portion of 
each panel. The mean profile (in arbitrary units) of 
the epoch is overplotted. The gray-scales plots are 22 
displayed in the radial velocity frame (the radial ve- 
locities are referred to the line laboratory rest wave- 
length) . 
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Fig. 5.— Photometry of WR 134 in 1992. Filled tri- 
angles: APT; open circles: SPM (Mexico) 1-channcl 
photometry; filled circles: SPM 2-channel photom- 
etry; open triangles: Maidanak Observatory; open 
stars: HIPPARCOS. The bottom panel is an en- 
largement of the insert in the central panel. The sizes 
of the symbols correspond to 2- a error bars. In order 
to guide the eye in the bottom frame, a sinusoid with 
a period of 2.27 days (ephemeris of McCandliss et al. 
1994) has been overplotted. 
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Fig. 6.— Photometry of WR 134 in 1992, folded with 
the V = 2.27 day period. Only the data plotted in 
the middle panel of Figure 5 have been used. The 
symbols are coded as in Figure 5. 
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Fig. 7.— Photometry of WR 134 in 1995. Filled 
triangles: APT; filled circles: SPM (Mexico); open 
triangles: Crimean Observational Station. The sizes 
of the symbols correspond to 2-cr error bars. 
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Fig. 8. — Auto-correlation matrices of He II A4686 for 
each epoch. Thick and thin contours indicate a nega- 
tive or positive correlation in the pattern of variability 
presented by the same line profile at different radial 
velocities, respectively. The lowest contour is drawn 
for a significant correlation at the 99.9 % confidence 
level. 
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Fig. 9. — Correlation matrices of He II A5412 with 
He II A4542, He II A4686, He II A4859, and C IV 
A5806 for epoch III. Thick and thin contours indi- 
cate a negative or positive correlation in the pattern 
of variability presented by two line profiles at differ- 
ent radial velocities (referred to the line laboratory 
rest wavelength), respectively. The lowest contour is 
drawn for a significant correlation at the 99.9 % con- 
fidence level. 
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Fig. 10. — Contours of constant ionization parameter 
log£ (in ergs cm s _1 ). The accreting, strongly ion- 
izing compact companion is located at (0, 17). The 
WR star is placed at the origin of this plot. 
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Fig. 11. — Upper and middle panels: predicted X-ray 
luminosities after allowance for wind absorption in 
the 0.2-2.4 and 0.2-4.0 keV bands as a function of the 
orbital phase <p, and for different values of the ioniza- 
tion parameter £: log£ = 2.1 ergs cm s _1 (solid line), 
l°g£ = 1-8 ergs cm s^ 1 (short-dashed line), and log£ 
= ergs cm s _1 (long-dashed line). The shaded aeras 
show the range of the ROSAT and Einstein count 
rates for WR 134, namely 0.46 ± 0.22 and 4.6 ± 1.6 
x 10 32 ergs s" 1 in the 0.2-2.4 and 0.2-4.0 keV bands, 
respectively (Pollock, Habcrl, & Corcoran 1995; Pol- 
lock 1987). The luminosities have been scaled to a 
distance of 2.1 kpc (van der Hucht et al. 1988). Note 
that the calculations do not consider the X-ray emis- 
sion intrinsic to the WR wind (generally ascribed to 
radiative instabilities), nor the scattering of the neu- 
tron star emission in the WR wind which may cause 
extra emission at X-ray eclipse. Lower panel: column 
density of the absorbing wind material in units of 10 23 
cm~ 2 , as a function of <p. 



Fig. 12. — Upper panel: solid line — the observed 
TVS; dotted line — the modeled TVS within the 
maximum-emission approach. Lower panel: solid line 
— mean epoch I profile of He II A4686; dotted line — 
the SEI model fit within the maximum-emission ap- 
proach; dashed line — the SEI model fit applying the 
minimum-emission approach. Note that the model is 
unable to fit the blucmost and redmost parts of the 
line profile because of partial blending and electron 
scattering effect, respectively. The projected velocity, 
w, is normalized to = 1900 km s _1 . 
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Fig. 13. — Sketch of the geometry adopted for the 
photoionization cavity. The azimuthal extension of 
the cavity, O, is in the direction perpendicular to the 
plane of the figure. 
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Fig. 14. — The observed (solid line) and modeled 
(dotted line) phase-dependent variations of He II 
A4686. All epoch I spectra are binned to 0.1 phase 
resolution; the bin mid-phase is shown in each panel. 




Fig. 15. — Variations of the P Cygni absorption 
troughs of He I A4471 and N V A4604 with phase. 
All epoch I observations are binned to 0.2 phase reso- 
lution: dashed- dotted line: <fi = 0.15-0.35; long-dashed 
line: <fi = 0.35-0.55; short-dashed line: <fi — 0.55-0.75; 
dotted line: <fi — 0.75-0.95; solid line: <j) — 0.95-1.15. 
Note that the intensity scale is different for the two 
panels. 
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Fig. 16. — The Stokes parameters Q and U of WR 
134 in 1985-1986 (Robert et al. 1989) plotted in 
phase according to the ephemeris of McCandliss et 
al. (1994). Left and right panels: data (with 2- a 
error bars) for the intervals HJD 2,446,362-2,446,371 
(1985) and HJD 2,446,704-2,446,735 (1986), respec- 
tively. 
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